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Abstract 

We have observed the Cygnus Loop from the northeast (NE) rim to the southwest 
(SW) rim using Suzaku in 10 pointings that is just north of previous XMM-Newton 
observations. The observation data obtained were divided into 45 rectangular regions 
where the width were configured so that each region holds 8000—12000 photons. The 
spectrum acquired from each region was fitted either with one-/cT e -component non- 
equilibrium ionization (NEI) model or with two-/cT e -component NEI model. The 
two-/cT e -component model yields significantly better fit in almost all the non-rim 
regions. Judging from abundances and flux, the high-/cT e -component (0.4— 0.8keV) 
must be the ejecta origin, while the low-/cT e -component (~0.3keV) comes from the 
swept-up matter. We found that swept-up matter shell is very thin in just southwest 
of center of the Loop. Together with previous observations, we estimate the diameter 
of this thin shell region to be 1°. We also found that the ejecta distributions were 
asymmteric to the geometric center: the ejecta of O, Ne and Mg were distributed 
more in the NE, while the ejecta of Si and Fe were distributed more in the SW of the 
Cygnus Loop. We calculated the masses for various metals and estimated the origin 
of the Cygnus Loop as the 12 — 15M core-collapse explosion. 

Key words: ISM: abundances - ISM: individual (Cygnus Loop) - ISM: supernova 
remnants - X-rays: ISM 
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1. Introduction 



The Cygnus Loop is a nearby (540pc: Blair et al. 2005) middle-aged supernova remnant 
(SNR). Its large apparent size (2°. 5 x 3°. 5: Levenson et al. 1997; Aschenbach & Leahy 1999) 
and high surface brightness enable us to study soft X-ray emission from the Cygnus Loop. 
The bright surface emission is mainly from interstellar medium (ISM) swept-up by forward 
shock wave. Miyata et al. (1994) observed the northeast (NE) rim of the Loop with Advanced 
Satellite for Cosmology and Astrophysics (ASCA) and revealed metal deficiency there. Since 
previous measurements indicate that the ISM around the Cygnus Loop has metal deficient 
abundances (Parker et al. 1967), they concluded that the plasma in the NE rim is dominated 
by the ISM. In the contrast of the rim region, ASCA detected Si, S, Fe rich plasma at the 
center portion of the Cygnus Loop (Miyata et al. 1998), which is thought to be the ejecta. The 
relative abundances of the ejecta support the idea that the Cygnus Loop was the result from a 
core-collapse supernova, and the progenitor mass is estimated to be 25M . 

Recent XMM-Newton observations (Tsunemi et al. 2007) across the Cygnus Loop from 
the NE rim to the southwest (SW) rim revealed ejecta distributed inside of ~0.85 R s of the 
Loop, where R s is the shock radius. The relative abundances inferred for the total ejecta are 
almost consistent with those expected for the core-collapse SN, whose progenitor mass is 13M 
(Katsuda & Tsunemi 2008). Suzaku observed just south of XMM-Newton observation path 
(Katsuda et al. 2008) and revealed asymmetric metal distributions of the Cygnus Loop. They 
found that Mg is distributed more in the NE, while Si, S, and Fe are distributed more in the 
SW of the Loop. In order to extent our knowledge on the ejecta distribution as well as the 
overlying ISM, we observed the Cygnus Loop from the NE rim to the SW rim using Suzaku 
observatory. Figure 1 shows our field of view (FOV) which cover just north regions of the 
XMM-Newton observation path. 

2. Observation and Data Screening 

The observations comprised ten points: NE3 , PI— P7, P9 and P10. We employed re- 
vision 2.0 of the cleaned event data since there was no sufficient flares in the light curve. For 
background subtraction, the spectra acquired from Lockman Hole were used. We selected the 
Lockman Hole data whose observation dates were close to those of our Cygnus Loop observa- 
tions. Since there was no photons above 3.0 keV after the background subtraction, the energy 
ranges of 0.3-3.0 keV and 0.4-3.0 keV were used for XIS1 (back-illuminated CCD; BI CCD) and 
XIS0, 3 (front-illuminated CCD; FI CCD), respectively (Koyama et al. 2007). We summarized 
Obs IDs, nominal points, observation dates, and effective exposure times after the screening in 
table 1. Figure 2 right shows a merged XIS1 three-color image. Red, green, and blue colors 
correspond to narrow energy bands of 0.52-0.70 keV (O VII Ka), 0.70-0. 85 keV (Fe L) and 
0.85-0.94 keV (Ne IX Ka), respectively. The scale R shows the distance from an observation 

2 



center in arcmin. We can see strong green color around the SW portion (i.e., P6, P7, P9), 
whereas red and blue colors are enhanced in the NE (i.e., PI— P4) regions and P10. 

3. Spatially resolved spectral analysis 

We divided the entire FOV into 2 parts, namely the NE part and the SW part. The NE 
part contains pointigs of NE3, PI, P2, P3, P4 and the SW part contains pointings of P5, P6, 
P7, P9, P10. The two parts were divided into total of 45 rectangular regions where their widths 
were configured so that each area holds 8000—12000 counts of photons in XIS0. The position 
angle of rectangular regions in the NE part and the SW part are set so that regions are parallel 
to each rim's shock front. The generation of the response matrix file (RMF) and the ancillary 
response file (ARF), are done by xisrmfgen (ver. May 2007) and xissimarf gen (ver. March 
2008). In order to reduce the effect of radiation damage of XIS, the data were taken by using 
the spaced row charge injection (SCI) method (Prigozhin et al. 2008) which recovers the energy 
resolution of radiation damaged CCD. The SCI method is not fully supported by CALDB at the 
writing phase of this paper, and our data clearly showed better energy resolution than that of 
the RMF generated by normal use of xisrmfgen. Therefore we searched the appropriate RMF 
by changing the date_obs parameter of xisrmfgen which represents the observation date. By 
changing the date_obs to earlier date, xisrmfgen can generate the RMF for earlier observation 
which has better energy resolution. We changed date_obs by every 3 months from September 
2005 to September 2007, which generate 9 different RMFs in total. Then the generated RMFs 
were used to fit the 10 sample spectra taken from red regions in figure 2. We applied model 
described in next section to test the RMFs in order to find the most appropriate RMF. Since 
the RMF where we set date_obs = December 2005 showed the least reduced-chi squared value 
in every spectrum, those RMFs were used to perform further spectral fit. 

3. 1 . One-kT e VNEI model 

The extracted spectra were first fitted by an absorbed non-equailibrium ionization (NEI) 
model with a single component [the wabs; Morrison & McCammon 1983 and the VNEI model 
(NEI ver. 2.0); e.g., Borkowski et al. 2001 in XSPEC vl2.4.01]. Free parameters are the 
hydrogen column density, 7V H ; electron temperature, /cT e ; the ionization time, r, where r is 
the electron density times the elapsed time after the shock heating; the emission measure EM 
(EM= / n e nn^, where n e and Tin are the number densities of electrons and hydrogens and dl 
is the plasma depth); the abundances of C, N, O, Ne, Mg, Si, Fe and Ni. We set abundance 
of S equal to that of Si and that of Ni equal to that of Fe. Abundances of the other elements 
were fixed to the solar values (Anders & Grevesse 1989). The black line in figure 3 shows the 
reduced-chi square value as a function of R. The one-component model gave us large reduced- 
chi square values in some regions; therefore, we added an extra component for our spectral 
fit. 
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3. 2. Two-kT e VNEI model 

We tried to apply this model by freeing all parameters, but since this did not give as 
reasonable values, we decided to fix certain values. The fcX e , r, and EM are free parameters in 
both components. is also a free parameter but shares same value between two components. 
Since abundances of the swept-up matter region is well known, we set the abundances of low- 
fcT e -component to those values determined in the NE region by Uchida et al. (2006). The 
abundances of O (=C=N), Ne, Mg, Si (=S), and Fe (=Ni) for the high-/cT e -component are 
set as free parameters. The fit improved significantly (> 99% based on F-test probability) in 
non-rim region (—70' < R < 65') by adding another component (see figure 3). On the other 
hand, the NE rim (R < 70') and the SW rim (R > 65') did not show notable improvement on the 
fit. Furthermore, parameters of the high-/cT e -component had large uncertainties in rim regions. 
Therefore, we employed two-/cT e -component model in the non-rim regions (— 70' < R < 65'), 
while one-/cT e -component model in the rim regions. Figure 4 shows spectra extracted from red 
regions shown in figure 2. The top 2 panels are fitted with one-component model (NE3 and 
P10) while the bottom 8 panels are fitted with two-component model. The best-fit parameters 
are summarized in table 2. The residuals between the data and the model in bottom 4 panels 
(P5, 6, 7, 9) of figure 4 show some structures around 1.2keV. We should keep in mind that 
the two-/cT e VNEI model still did not give us acceptable fits from the statistical point of view, 
which suggests that our model is too simple. 

Figure 5 shows fcT e , r and metal abundance as a function of R. We found that the 
high-/cT e -component shows distinctly higher temperature and metal abundance compared to 
low-/cT e -component. The temperature of the low-/cT e -component stays almost constant (~ 0.3 
keV) while the temperature of the high-/cT e -component has a peak value of 0.8 keV around 
R = —60' and decrease towards the SW. We also found asymmetry in metal abundance of 
high-/cT e -component. The abundance of 0[=C=N], Ne and Mg are higher in the NE part 
(R < 0') while Si is in the SW (R > 0'). The flux of each component is shown in top left 
panel of figure 5. The flux of the low-/cT e -component is dominant in the NE part and it 
also shows shell brightening in both NE and SW rims, therefore, we confirm that low-kT e - 
component surrounds the Loop and the high-/cT e -component fill its interior . All of these facts 
support our assumption where low-/cT e -component originates from the swept-up matter and 
the high-/cT e -component originates from the eject a. 

4. Discussion and Conclusion 

We observed the Cygnus Loop from the NE to the SW with Suzaku in ten pointings. 
Our FOV covers more to the north of previous observation carried out by XMM-Newton and 
Suzaku. Dividing the entire FOV into 45 rectangular regions, we extracted spectra from all the 
regions, and performed spectral analysis to them. For spectral fit, we employed models used 
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by the previous observation (Katsuda et al. 2008; Tsunemi et al. 2007) and obtained similar 
results. The one-component VNEI model showed fairly good fit in the rim of the Loop but not 
in the non-rim part of the Loop. The fit of non-rim part improved significantly by applying 
two-component VNEI model. Judging from the abundanes and the flux distribution of each 
component, the high-/cT e component must be the ejecta, while the \ow-kT e component comes 
from the swept-up matter. 

4-1. Swept-up Matter Distribution 

The temperature for the swept-up matter component is significantly lower than that for 
the ejecta component. This temperature is similar to that obtained for the rim of the Loop, 
where we expect no contamination of the ejecta (e.g., Miyata et al. 2007). Therefore, we 
believe that we surely separated the X-ray emission of the ejecta inside the Loop from that of 
the surrounding matter. The EM distribution of the swept-up matter is inhomogeneous in our 
FOV, as shown in top right panel of figure 6. The shell of the swept-up matter seems to be thin 
around 10' < R < 60' relative to that in the NE part (R < 0'). The flux of swept-up matter also 
shows this tendency; its flux in the SW part is about a third of the NE part. Such trend is also 
reported by XMM-Newton (Tsunemi et al. 2007) and previous Suzaku observations (Katsuda et 
al. 2008). Tsunemi et al. (2007) divided XMM-Newton observation into north path and south 
path and found that this thin shell region to be 5' in south path and 20 / in north path. They 
estimated this thin shell region to have diameter of 1° and centering (20M9 m ll s , 31°50 / 20 // ). 
Our observation path goes though right in middle of this region and our result was consistent 
with their prediction. Since the EM of this region is also about \ — | of the NE part, and the 
EM of ejecta is stronger compared to the NE part, there might be a blowout in the direction 
of our line of sight, just like a blowout in the south of the Loop (Uchida et al. 2008) and we 
are only seeing the one side of the shell. We detected no ejecta emission from the NE rim 
[R < —70') and the SW rim (i? > 65') of the Loop, which is about 15% of our FOV. The ejecta 
occupy the major part of the inner side of the Loop while the swept-up matter surrounds it. 
This structure is identical to those in previous results (Miyata et al. 1998; Katsuda & Tsunemi 
2008). 

4.2. Metal distribution in Ejecta 

We have calculated the EM of various heavy elements in the ejecta such as 0[=C=N], 
Ne, Mg, Si, and Fe [=Ni]. Figure 6 shows the EM distribution for these elements as a function 
of R. The red marks show the results from our FOV while the green and the blue marks show 
the results from XMM-Newton north path and south path. The comparison of these results 
show some similarity near the both rim of the Loop (R < — 60', 40' < i?), but shows discrepancy 
in the center portion of the Loop. This is especially notable in Si and Fe. Since our FOV is 
further north from the geometric center of the Loop compared to XMM-Newton observation, 
these discrepancies suggests onion layer structures at a SN time of explosion. For Si and Fe, 
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some structures can be seen in both Suzaku and XMM-Newton observations, such as the bump 
around R = -50 / and the drop at / < R < 20'. 

From the result from our FOV, the NE part and the SW part have different characteristic 
of metal distributions. We calculated the mass ratio in the NE part and those in the SW part 
from Figure 6. They are 0~ 4.7, Ne^ 3.2, Mg~ 2.3, Si~ 0.36, Fe~ 0.62 assuming unity filling 
factor. The O-Ne-Mg group is heavily distributed in the NE part of the Loop by factor of 
2 — 5. In contrast to the O-Ne-Mg group, the Si-Fe group is distributed more in the SW part 
of the Loop by factor of 2 — 3. These asymmetries can also be seen in all of the abundance 
distribution (figure 5) except for Fe. However EM of the elements is the clearer representation of 
the amounts of elements, therefore we believed that Fe is also distributed asymmetrically more 
to the SW. This asymmetry is quantitatively consistent from previous observation (Katsuda 
et al. 2008). A natural explanation for the asymmetry is an asymmetry at the time of the 
SN explosion of the Cygnus Loop. Recent theoretical models describe asymmetric supernova 
explosion resulting from hydrodynamic instability (Burrows et al. 2007). 

We calculated the mass of ejecta to be 21M from EM. We estimated the plasma depth 
to be 26pc, unity filling factor and n e = 1.2n#, although the fossil ejecta might be deficient 
in hydrogen. If this is the case, the total mass of the ejecta reduces to ^12M However this 
value strongly depends on the plasma structure of the ejecta. Therefore we used SN explosion 
model to calculate the mass. In order to compare our data with the SN explosion models, we 
calculated the ratios of elements relative to O. Figure 7 shows the number ratios of Ne, Mg, 
Si and Fe relative to O of the ejecta component. We plotted the core-collapse model (Woosley 
& Weaver 1995) for the various progenitor mass and Type la supernova model (Iwamoto et al. 
1999) for comparison. Since the distribution characteristic of elements showed large difference 
between the NE and the SW, we also plotted the number ratios calculated from each part. The 
orange line represents the number ratios calculated from all of observations including XMM- 
Newton and previous Suzaku observations (Tsunemi et al. 2007, Katsuda et al. 2008). The 
number ratios from the SW part shows fairly good agreement with Type la models in Mg and 
Si but that the number ratios of Ne and Fe show large inconsistency from the model. On the 
contrary, the number ratios of Ne, Mg, Si, from NE are in agreement with 15M model, but 
Fe is about 4 times higher in our data. Because of the asymmetric distributions of ejecta, the 
number ratios of heavy elements give very different values by choosing which part of the Loop 
we use. The number ratios calculated from our FOV and all FOV (2 Suzaku and 1 XMM- 
Newton) are both in good agreement with that of 12M model. However, just like any other 
data, the number ratio of Fe shows disagreement. It is possible that the Fe is overabundant 
because all FOV are across the center of the Loop, where we expect to be Fe rich region. 
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Fig. 1. ROSAT HRI image of the entire Cygnus Loop. The Suzaku FOV (NE3, PI, P2, P3, P4, P5, P6, 
P7, P9, P10) are shown as white rectangles. Dotted circles and rectangles represent previous XMM- Newton 
and Suzaku observations. 



Table 1. Information of observations of the Cygnus Loop and Lockman Hole 



Obs. ID Coordinate (RA, DEC) Obs. Date Effective Exposure 



Cygnus Loop 



500022010 (NE3) 


313.746, 32.188 


2005.11.29 


12.2 ks 


501012010 (PI) 


313.510, 31.975 


2007.11.13 


9.8 ks 


501013010 (P2) 


313.265, 31.779 


2007.11.14 


16.4 ks 


501014010 (P3) 


313.032, 31.574 


2007.11.14 


7.5 ks 


501015010 (P4) 


312.799, 31.369 


2007.11.14 


18.3 ks 


501016010 (P5) 


312.547, 31.180 


2007.11.15 


19.3 ks 


501017010 (P6) 


312.297, 30.991 


2007.11.11 


28.7ks 


501018010 (P7) 


312.078, 30.776 


2007.11.12 


21.0 ks 


501019010 (P9) 


311.809, 30.603 


2007.11.12 


16.2 ks 


501020010 (P10) 


311.566, 30.407 


2007.11.13 


14.6 ks 




Lockman Hole 






100046010 (for NE3) 


163.4063, 57.6108 


2005.11.14 


49.3 ks 


102018010 


162.9257, 57.2581 


2007.05.03 


68.9 ks 




Fig. 2. Left: Surface brightness map of our FOV. Right: Three-color image of ten XIS FOV (Red: O VII 
Kct, Green: Fe L, Blue: Ne IX Ka). The data were binned by 8 pixels and smoothed by a Gaussian kernel 
of (7 = 25" The scale R shows the distance from observation center in arcmin. The effects of exposure, 
vignetting, and contamination are corrected for both figure. 



9 



Table 2. Best-fit parameter for example spectra. NE and SW (P10) are fitted with one-component VNEI model while the others are fitted with two-component VNEI 
model 



Elements NE3 PI P2 P3 P4 P5 P6 P7 

Free abundance component 

n H [xl0 20 cm- 2 ] 3.8±0.7 4.1 ±0.8 3.3±0.5 4.5±0.5 3.1±0.6 3.2±0.5 3.2±0.5 3.6±0.4 

/cT e [keV] 0.28 t° ° l 0.74 ±g;g? 0.66 ±0.01 0.63±0.02 0.55±0.01 0.53±0.02 0.48± 0.01 0.41 ±0.01 

C 0.3 ±0.1 =0 =0 =0 =0 =0 =0 =0 

N 0-13±o.o5 =0 =0 =0 =0 =0 =0 =0 

O O-lOto.04 1-10 ±0.07 0.73±0.06 1.00 ± 0.07 0.52 ±0.04 0.15 ±0.02 0.11 ±0.01 < 0.1 

Ne 0,20±g;g| 0.85 ±0.08 0.76±0.07 0.71 ±0.07 0.18±0.03 0.15±0.03 0.14 ±0.03 < 0.1 

Mg 0.12 ±S:8? 0.51 ±0.07 0.47±0.06 0.34 ±0.07 0.13 ±0.03 0.12±0.02 0.16 ±0.03 0.12 ± 0.03 

Si 0.15 ±0.09 0.70 ±0.1 0.8±0.1 1.5 ±0.2 1.1 ± 0.1 1.7±0.1 2.6±0.2 1.3±0.1 

Fe(=Ni) 0.14±g;J 4 0.53 ±0.05 0.66 ±0.04 0.82 ±0.04 0.47±0.02 0.66±0.02 0.88 ±0.02 0.53±0.01 

r[xlO n cm- 3 s] 0.99 ±0.04 0.83 ± 0.06 1.00 ±0.04 1.3 ±0.1 1.92±g;J 9 1.53 ±g;?| 1.33 ±g;2 7 2.40 ±g;| 

EM[xl0 19 cm- 5 ] 4.20±g;| 0.05 ±0.01 0.046 ±0.009 0.03±g;gJ 9 0.10±g;g? 3 0.077 ±0.003 0.076±g;gg| 0.178±g;gJ 6 

Fixed abundance component 

kT e [keV] 0.27± 0.01 0.26 ±0.01 0.26±0.01 0.27±0.01 0.34 ±0.01 0.28 ±0.01 0.27±0.01 

Abundances (fixed to those determined for the NE rim of the Cygnus Loop)"'" 

r[xlO n cm- 3 s] . 0.56 ±0.04 0.97 ±0.04 1.02 ±0.07 0.72 tg;^ 0.54 ±g;g§ 0.92 ±g:g| 1.49±0.08 

EM[xl0 19 cm- 5 ] 1.6 ±0.2 1.4 ±0.2 1.1 ±0.1 0.787t°;f - 17 ^o.oo7 °- 20 ^o.oi - 21 ^o.oi 

* Other elements are fixed to solar values. 

The values of abundances are multiple of solar value. 

The errors are in the range A% 2 <2.7 on one parameter 

tC=0.27, N=0.11, 0=0.11, Ne=0.21, Mg=0.17, Si=0.34, S=0.17, and Fe (=Ni)=0.20 (Uchida et al. 2006) 



P9 



SW(P10) 



3.5 ±0.5 
0.40 ±0.01 
=0 
=0 
0.11 ±0.02 
0.12 ±0.02 
0.10 ±0.03 

1.5 ±0.2 
0.56 ±0.02 



3.2 ±0.6 
0.32 ±0.01 
0.38±0.08 

u - iO -0.05 

16 +0 - 08 

u - iD -0.06 



0.21 



+0.1 



3.3 



+0.8 
-0.9 



0.08 

0.14 ± 0.02 
0.31±0.1 
0.19±0.1 
1.6 ±0.05 



0.11 ±0.01 0.37 



+0.3 
-0.2 



0.33±0.01 



0.71 ±0.04 
n iQ+o.oi 

U ' i ^-0.02 
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0.5 1 2 0.5 1 2 

Energy (keV) Energy (keV) 



Fig. 4. X-ray spectra extracted from the red regions in Figure 2. The best-fit curves of two-component 
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Fig. 5. The fcT e , r values and 
i?. The black marks show 



abundance of elements as 
the parameters of low-kT e 



a function of 
for comparison. 
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-50 50 o -so 50 




Fig. 6. The top left panel shows the EM (J n e riHdl) in 10 19 cm -5 . The other panels are EM of elements 
(J n e nxdl) in 10 14 cm -5 . The results from our Suzaku FOV is shown in the red mark while the green and 
the blue mark shows the results from XMM-Newton north path and south path (Tsunemi et al. 2007). 
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Ne Mg Si Fe 



Fig. 7. Number ratios of Ne, Mg, Si, Fe, relative to of the high-/cT e -component for the entire FOV 
(black line). Black dotted line is a number ratio for NE part while black dot-dash line is for SW part 
of the Loop. The orange line is a number ratio for all of the 2 Suzaku observations and XMM- Newton 
observations (Katsuda et al. 2008; Tsunemi et al. 2007). Dotted red lines represent the CDD1 and W7 
Type la SN models of Iwamoto et al. (1999). Dotted blue, light blue, magenta, and green lines represent 
core-collapse models with progenitor masses of 12, 13, 15, 20 M respectively (Woosley & Weaver 1995). 
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